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The NPS computer simulation model was used to investi- 
gate the sputtering by 1.0 keV argon bombardment of clean 
and oxygen reacted monocrystalline titanium, vanadium and 
niobiiim. Variations in yield, energy and angular distri- 
butions and in particular, multimer formation, were studied 
as a function of oxygen density and location. Simulation 
results show a significant decline in substrate yield as 
the oxygen coverage is increased, regardless of location. 
Fxirther, there is a marked preference for multimer forma- 
tion by lattice fragmentation, rather than by recombination, 
for all three metals. The percentages of multimers formed 
by fragmentation were found to increase with increasing oxy- 
gen density for all three metals. Previous results on face- 
centered cubic crystals concerning the yield per layer and 
the effects of channeling were confirmed for body-centered 
cubic and hexagonal closed-packed crystals. For compari- 
son, simulation runs were repeated at ion energies of 0.5 
and 2.0 keV for a titanium with oxygen atoms located in 
three-fold sites with a C(2x2) density. 
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I . BACKGROUND 



A. OVERVIEW 

Sputtering is the ejection of atoms from a target by 
bombardment with an energetic projectile ion. The incoming 
ion collides with atoms of the solid, thereby losing a 
portion of its energy and moment\am. If the energy gained 
by a struck atom exceeds the binding energy holding it in 
position in the lattice a primary recoil atom will be 
created. This recoil atom will then collide with other atoms 
in the material, distributing its energy via a collision 
cascade. An atom at or near the surface will be sputtered 
if the energy associated with the normal component of its 
velocity exceeds the surface binding energy of the material. 

The phenomenon of sputtering can be subdivided into two 
principle categories. Transmission sputtering is the ejec- 
tion of atoms from the rear of a thin target following pene- 
tration by the ion. Sufficient energy must be transported 
through the target to allow atoms to overcome their binding 
energy and escape. Back-sputtering is the ejection of atoms 
from the front of the target and depends upon sufficient 
energy being deposited in the surface layers to allow 
ejection. This thesis will be concerned with the latter. 

First discovered by Grove in 1853 and Faraday in 1854 
as the deposition of metal atoms on the glass walls of a 
gas discharge tube, it was not until 1902 that Goldstein 
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[Ref. 1] presented evidence that the sputtering effect was 
caused by positive ions of the discharge impacting on the 
metal cathode. Since its principle effect was to contaminate 
experimental environments and to erode apparatus, the majority 
of early work centered around ways to control or eliminate 
what was considered to be an undesirable side effect. The 
necessity of increasing the life expectancies of high voltage 
vacuum tubes and the recognition of the possible applications 
of sputtering to thin film coatings gradually spurred inter- 
est in understanding the mechanics involved. 

Early theoretical work was hampered by a paucity of 
experimental data. Reproducibility of results was also a 
major problem due to an incomplete understanding of the 
factors involved. Since the sputtering yield, Y (defined as 
the number of atoms ejected per incident ion) , is critically 
sensitive to ambient pressure and surface contamination, as 
well as the flux density and angle of incidence of the beam, 
progress was dependent on the development and refinement of 
experimental techniques . 

B. EARLY EXPERIMENTAL WORK 

The effect of pressure on yield was first demonstrated 

by Penning and Moubis [Ref. 2] . They found that an increase 

in pressure increased the collision frequency between the 

ejected atoms and surrounding gases, resulting in the back- 

scattering of the escaping atoms back onto the target 

-5 

surface. By keeping background pressure below 10 Torr, 
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they were able to obtain reproducible results for ion ener- 
gies in excess of 500 eV. Below 500 eV, the scouring action 
of the incident beam was insufficient to maintain a clean 
surface. 

The elimination of adsorbed gases and surface oxide 
layers as a critical factor in determining the true yield 
was described by Arifov, et al., in 1963 [Ref. 3]. Yonts 
and Harrison also presented evidence [Ref. 4] that surface 
recontamination from background gases was a significant 
factor in quantitative sputtering yield measurements. 
Criteria were therefore developed for conducting 'clean' 
experiments . 

In addition to pressure and surface cleanliness, the 
following conditions must be met in order to obtain repro- 
ducible yields: 

(1) The ions must strike the target at a known angle 
of incidence. 

C2) The energy spread of the beam must be small. 

(3) All ions must be uniformly charged and mass separated 

(4) On monocrystalline surfaces the lattice orientation 
must be adequately described. 

C. THEORETICAL DEVELOPMENT 

Pioneering work on the concept of an individual sputter- 
ing event on an atomic scale was conducted by Stark [Ref. 5] 
Applying the conservation laws governing elastic collisions, 
he proposed two theories to explain known experimental 
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results. One, the "Hot-Spot” model, considered the sputtered 
atoms to be the result of evaporation of target material 
from an atomically small region due to ion induced local 
heating. In his "Collision" model, sputtering events were 
the result of a series of binary collisions initiated by 
a single ion. His models were modestly successful in pre- 
dicting the energy distribution of the sputtered particles. 
Further refinements on the hot-spot model were attemp'ted by 
von Hippie [Ref. 6], but only after the demonstration by 
Wehner [Ref. 7] of the effects of crystal structure on the 
yield did it become apparent that local heating alone could 
not explain the mechanisms involved in sputtering. 

Wehner 's observation of spot patterns (angular distri- 
bution patterns characteristic of the crystallographic proper- 
ties of the substrate) obtained from sputtering of mono- 
crystalline targets revived interest in the collision approach. 
In 1952 Keywell [Ref. 8] made a first attempt to formulate 
Stark's collision model in terms of the neutron transport 
model developed for nuclear work. Harrison then developed 
a theory based on probability concepts as expressed by the 
idea of collision cross-sections [Ref. 9] . However, due to 
the large number of unknown parameters, these theories were 
limited in their quantitative application. 

As a further attempt to account for the angular distri- 
bution of ejecta, Silsbee [Ref. 10] proposed a focused 
collision model which allowed the transport of momentum in 
crystals along preferred directions. Available experimental 
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results showed that the yield from monocrystals depended 
sensitively on the crystallographic orientation of the inci- 
dent beam. This was explained in terms of the "holes" seen 
in crystal models when viewed from different planes and 
led to the development of various "Transparency" models. 

While focusing has been shown to be a contributing factor 
in cascade development within a crystal , it is inadequate 
to explain Wehner spots. 

As an alternative to the focusing collision theory, 
Lehmann and Sigmund [Ref. 10] proposed a model based on the 
Boltzman transport equations that required the target surface 
to have an ordered structure, but not long straight rows of 
atoms intersecting the surface as in the focusing model. 
Emphasis is placed on the role of the surface structure and 
in particular, surface binding energies in monocrystalline 
sputtering. As an ion strikes the surface a collision cas- 
cade is created via binary interactions. An atom in the 
cascade will sputter if its momentum component normal to 
the surface has sufficient energy associated with it to 
overcome the surface potential barrier. Thompson [Ref. 12] 
expanded on this idea and proposed that the surface attrac- 
tion for the escaping particle causes a refraction of its 
velocity vector away from the noirmal, resulting in a dis- 
torted angular distribution of ejected particles. Signifi- 
cant agreement with experimental results was found when 
projectile energies were much greater than the surface bind- 
ing energy. This approach also provided correct information 
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concerning both the angular and the energy distributions 
from polycrystalline targets. Thompson's [Ref. 12] work 
also incorporated the idea of focusons (momentum transport 
without mass displacement) into the theory. Robinson (Ref. 

I 

13] has taken a similar approach while accounting for 
channeling. 

While these theories have provided mathematical tools 
to aid experimentalists, they have been rather unsuccessful 
in predicting yields at low projectile energies and from 
single crystal targets. Binary collision models have thus 
far been unable to account for the predominance of sputtered 
atoms from surface layers of the target. 

In general, these theories have attempted to reduce what 
is essentially a many-body, multiple interaction process 
into forms that are analytically tractable. An assortment 
of refinements followed, attempting to explain observed 
behavior. Full scale attack on the multiple interaction 
problem had to await the arrival of high speed computers with 
sufficient storage to handle the hundreds or thousands of 
atoms in a cascade. 

D. COMPUTER SIMULATION MODELS 

Paralleling the theoretical and experimental progress in 
the understanding of sputtering was the effort to model the 
basic mechanisms through the use of computer simulation. In 
1960 Gibson, Goland, Milgram and Vineyard [Ref. 14J developed 
a computer model to simulate the motion of atoms in a copper 
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monocrystal following tho impact of an ion on a single 
ts-rget atom. Atomic interactions were treated as binary 
collisions and the resulting cascades were developed by 
applying Newton's equations of motion. Ion penetration 
studies by Robinson and Oen [Ref. 15] using a similar 
program led to the discovery of channeling which was later 
experimentally verified. 

In 1967, based on earlier preliminary studies, Harrison, 
Levy, Johnson and Effron [Ref. 16] used a computer to simu- 
late the bombardment of a single copper crystal by argon. 

An important finding of this work was that the yield was 
almost entirely confined to the first three layers of the 
top surface. Subsequent refinements of the model [Refs. 17- 
20] included the addition of potential functions that allowed 
for attractive, as well as repulsive interactions, creating 
a dynamically stable crystal. Harrison's results showed that 
focusons were not a major contributor to the yield. 

In 1978 Garrison, Winograd and Harrison [Ref. 21J pub- 
lished the results of a comprehensive study of atomic and 
molecular ejection from a copper crystal with adsorbed oxy- 
gen atoms . Their simulation of the system provided a de- 
tailed picture of the mechanisms of ejection and enabled 
them to determine whether molecules were ejected as a cluster 
from the surface or combined together in the space above the 
target. They were also able to determine the effects of 
adatom placement on molecule formation. Such results, when 
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coupled with laboratory data, should lead to a better under- 
standing of actual adatom locations on chemically reacted 
surfaces . 

While much work remains to be done, particularly with 
respect to constructing interatomic potentials that accurately 
represent the physical processes involved, computer simula- 
tion permits the study of sputtering phenomena with a preci- 
sion and flexibility that is impossible to duplicate in the 
laboratory. While it cannot replace actual experimentation, 
it can lead to a better understanding of experimental re- 
sults and provide a testing ground for theoretical concepts . 
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II. OBJECTIVES 



The primary objective of this thesis is to study the 
sputtering of titanium, vanadium and niobium monocrystals 
using the NFS computer simulation mpdel, QRAD. In addition, 
changes in the yield and multimer formation are studied when 
the surface is reacted with oxygen. 

Yu [Refs. 22,23] has conducted comprehensive research 
on clean and oxygen reacted titanium, niobium and vanadium 
in an attempt to link several earlier theories on the 
mechanisms responsible for ionic multimer formation. In 
the Lattice Fragmentation Model first proposed by Benning- 
hoved [Ref. 24] molecular species Cmultimers) are emitted as 
complete fragments of the surface lattice during ion bombard- 
ment. In the Recombination Model [Refs. 18,25] only single 
atoms are ejected. Molecules are then formed through the 
attractive interaction between sputtered atoms from the 
same sputtering event, establishing their identities above 
the surface. The charge state of the molecules is not taken 
into account in these models. 

Yu has reported experimental evidence [Ref. 23] that the 
emission mechanism for a sputtered molecular species may be 
a function of its charge state. He proposed that for low 
oxygen coverages (< one monolayer) the formation of positively 
charged titanium and niobiiam molecules (MO^, MO^) favored 
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the re coinbi nation model, while the formation of negatively 
charged species (MO , ^^2 ^^ 3 ^ followed the fragmenta- 

tion process in all three metals. In the case of vanadium 
it appeared that the positive oxide species were not formed 
exclusively by recombination, but by a combination of the 
two methods . 

Formation of positively charged ions was also found to 
be strongly dependent on oxygen coverage. When the niomber 
of oxygen atoms neighboring a metal atom increases , the 
frequency of sputtering an and 0^ adjacent to each other 
naturally increases . Consequently the percentage of MO^ 
produced by fragmentation will increase. 

The formation of negative molecular ions by lattice 
fragmentation was found to be insensitive to oxygen coverage 
due to the kinetics involved in the formation proces. A 
molecule formed by recombination is usually left with excess 
energy. If this energy approaches the electron affinity 
of the molecule the negative ion will attempt to stabilize 
itself by ejecting an electron. If the excess energy is 
much greater than this, one approaches the dissociation 
energy, and the molecule will be prevented from forming at 
all. Thus the formation of a stable negative ion is most 
favorable when the components are close neighbors on the 
surface and eject with small relative velocities as in the 
lattice fragmentation process. Increasing the oxygen cover- 
age will only increase the probability of forming higher 
order oxide multimers . 
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Some discrepancies were noted in Yu's results for posi- 
tive ions and were attributed to the recombination of metal 
atoms and oxygen from adjacent sites. Thus the results were 
indistinguishable from fragmentation at higher oxygen 
coverages. All nine of the negative molecular ion species 
investigated by Yu showed strong evidence of formation by 
fragmentation. Since the energy balance for recombined 
negatively charged ions indicates that they tend to neutralize 
by electron ejection, this leaves fragmentation as the 
principal soxirce of molecule formation. 

It should be noted that Yu's experimental procedure 
utilized an ion beam with a 45 degree angle of incidence to 
the polycrystalline surface. The mass spectrometer detector 
was also positioned at 45 degrees with respect to the surface, 
at right angles to the ion beam. Thus only a portion of the 
total yield was sampled. The data from the simulation studies 
that follow are based on the total yield from ion bombardment 
perpendicular to monocrystal surface. These data include 
neutral atoms, as well as ions, which are not collected by 
the mass spectrometer. 

While determination of the charge state of molecules is 
not possible with the computer model, the various mechanisms 
of ejection can be studied in detail to determine whether 
recombination or fragmentation is the major contributor to 
the multimer yield. Ejection times of the various constitu- 
ents are also available for correlation. 
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Additionally the energy and angular distributions of 
ejected particles, which can be related to experimental SIMS 
data, will be studied. Since the actual sites of adatom 
adsorption have yet to be determined, the effects of oxygen 
coverage and location on yield and multimer formation will 
also be investigated. 
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III. MODEL DEVELOPMENT AND SIMULATION 



A. COMPUTER MODEL 

9 

QRAD is a computer simulation program designed to track 
the motions of atoms in a microcrystal following ion bombard- 
ment. Trajectories are calculated by solving Hamilton's 
equations of motion using an average force method [Ref. 26] . 

It is a multiple interaction (MI) simulation following time 
step logic. Once the atomic masses, crystal characteristics, 
interatomic potentials and ion energy have been initialized, 
the program will develop and follow the resulting collision 
cascade until the energy of the atoms reaches a predetermined 
minimum. (See reference 24 for a comprehensive discussion 
of sputtering simulations.) 

Initially an ion of given mass is projected toward a 
surface impact point with a specified energy and angle of 
incidence. As the cascade develops, the positions and veloci- 
ties of the ion and lattice atoms are calculated at the end 
of each timestep. A timestep is determined by the velocity 
of the most energetic particle so that its displacement does 
not exceed 0 J. lattice unit CLU) To expedite calculations, 
particles are classified as either moving or non-moving de- 
pending on their energies and the forces exerted on them 

^The lattice unit is a convenient method of expressing 
distances when dealing with different crystal types . An 
LU is equal to the lattice constant 'a' divided by root two 
for hep crystals and is equal to a/2 for bcc crystals. 
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by noighbor atoms . Atoms which rise above the surface with 
sufficient momentum are reclassified and assigned to a 
tentative list of ejected atoms. A trajectory is terminated 
when the energy of the most energetic atom is such that the 
probability that either it, or other atoms, will escape the 
surface is very small. Final positions and velocities are 
then recorded. Multiple trajectories, simulating bombard- 
ment by an ion beam, can be run by computing trajectories 
for a set of impact points chosen to sample a representative 
area of the crystal surface as a function of its symmetry. 
Data on the sputtered atoms is then analyzed by separate 
programs (AN83 and PLOTAN) to determine multimer formation, 
ejection angle information, etc. 

As mentioned previously, multimers may result either 
from lattice fragmentation or by recombination above the 
surface. A multimer is considered formed when its constitu- 
ent atoms are in a favorable spatial position with suffi- 
ciently low relative kinetic energy to permit bonding. The 
total energy of the resulting system must be negative and 
is an indication of the group's binding stability. Absolute 
multimer yield is strongly dependent on the pair-potential 
in use and the range specified for atomic interactions 
(R^) , however the mechanisms for cluster formation have been 
found to be insensitive to the choice of the potential 
function [Ref. 27] . 

Global results, suitable for comparison with experimental 
data, are obtained by averaging the yields from each run 
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over a number of trajectories. While the number of trajec- 
tories is quite small (100-200, depending on symmetry 
considerations) , the results have been found to vary by 
less than 10% when compared with other independently chosen 
sets of impact points [Ref. 28] . These global results are 
used in determining the number of atoms ejected per single 
ion (ASI yield) , energy and angular distributions, as well 
as the spatial distribution of yield from the impact zone. 

The remainder of this chapter is devoted to the develop- 
ment of the potential functions and other parameters used 
by QRAD. Adatom nomenclature and placement are discussed 
in Section D. Trajectory sets and crystal impact points are 
covered in Sections F and G. The chemical properties of 
titanium, vanadium and niobixam are covered in the following 
sections . 

B. TITANIUM 

Titanium is a dimorphic transition metal ranking ninth 
in abundance of the elements making up the earth's crust. 
Below about 880 ®C it exists as a silvery gray hexagonal 
close-packed (hep) crystalline material, while above 880 °C 
it undergoes a phase transformation to a face— centered 
cubic (fee) structure. The important physical and chemical 
properties of titanium are listed in Table 1. 

Titanixjm reacts with oxygen to produce four well defined 
oxides, the monoxide TiO, the sesquioxide Ti202/ the dioxide 
or titanic acid Ti 02 and the trioxide or pertitanic acid 
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TiOs* Ti 20 g, Ti^O^ and Ti^ 0 j ^2 have also been reported, but 
may in fact be metastable combinations of the others [Ref. 
29] . 

At least two distinct forms of oxygen coverage have been 
observed on Ti(OOl) surfaces. 

1. Type 1 begins to form at low oxygen exposures of 

2 

approximately 1.0 L and is characterized by a low work 
function, very low electron stimulated desorption (ESD) 
and photon stimulated desorption (PSD) yields and a surface 
binding energy of 5.9 eV as determined by UV photoemission 
spectroscopy (UPS). Plateau, et al., had previously deter- 
mined this value to be 5.7 ± 0.2 eV, also using UPS and 
ESCA [Ref. 38]. 

2. As the oxygen exposure is increased (~50 L) a second 
type of oxygen coverage is observed which has a higher work 
function, high PSD yield and additional binding energy in 
UPS. The effects of oxygen adsorption appear to saturate 
at about 70 L. 

Hanson [Ref. 37] concluded that type 1 may correspond 
to oxygen adatoms lying below the surface plane or in high 
coordination sites, possibly the three-fold hollows, which 
would account for the low desorption rate. 

Jonker [Ref. 39] has also observed two states of oxygen 
on the surface. The first is detected at low coverages and 
produces a well-ordered P(.2x2) low energy diffraction 

9 — 6 “*4 

1.0 L (Langmuir) = 10 Torr-sec = 1.33x10 Pa-sec. 
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(LEED) pattern. At higher oxygen exposures a second state 
emerges which has a disordered appearance. This was ex- 
plained by assuming that the sticking probability coefficient 
for oxygen is relatively constant up to a quarter monolayer 
coverage. At this point interference from previously 
adsorbed oxygens would cause it to change abruptly. The 
sticking coefficient would then remain constant for increas- 
ing exposures until the density of oxygen favored a change 
in the chemical bonding, resulting in the formation of a 
more stable metal oxide in the upper layers of the crystal. 

Kawasaki, et al. [Ref. 40], have also reached the conclu- 
sion that oxygen is first adsorbed on the surface at low 
doses, but as the exposure is increased, it is slowly ad- 
sorbed into the surface. This corresponds to the formation 
of an oxide, and is accompanied by changes in the crystal 
structure. Their results also indicate that at the lower 
exposures oxygen occupies either a doubly bonded A-^top or 
a two-fold bridge position. The A- top position is also sup- 
ported by PSD experiments [Ref. 41] in which the oxygen 
atoms were observed to desorb normal to the surface, imply- 
ing that the Ti-0 bonds are perpendicular to the surface. 

Experimental work by Singh, et al. [Ref. 43], confirms 
the formation of an oxide as saturation is reached, with 
oxygen atoms incorporated extensively into the metal lattice. 
(Here saturation is defined as the point where increasing 
the oxygen exposure no longer has an effect on the metal 
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crystal, i.e., the sticking probability is equal to zero.) 

Since the sputtering yield comes mainly from the top three 
layers, modeling of an oxygen saturated surface would have 
to be accomplished by using a titanium oxide crystal rather 
than titanium with oxygen adatoms. 

In view of the uncertainty of the actual locations of the 
oxygen atoms, a complete set of trajectories was run for the 
A-top, two-fold (bridge) and three-fold locations, in addi- 
tion to those on the clean surface. Figures 1 through 8 
present the various adatom locations on the metal crystal 
surface. (See Section D for a description of adatom notation.) 

C. VANADIUM AND NIOBIUM 

Vanadium and niobium are also transition metals, but 
unlike titanium, they have a body-centered cubic (bcc) 
structure. Their (110) planes, while not as close packed, 
do however have a hexagonal arrangement of atoms similar to 
the titanium (001) plane (Fig. 9) . Their physical and 
chemical properties (Tables 2 and 3) are also quite similar 
to titanium, hence their sputtering mechanisms should bear 
some correlation to that of titanium. 

Haas, Jackson and Hooker [Ref. 44 j in their DEED experi- 
ments on niobixam, tantalum and vanadium have found that oxygen 
adsorption is the same on all three metals and is a function 
of oxygen exposure. At low oxygen coverages the oxygen 
patterns clearly indicate a six-fold symmetry characteris- 
tic of the substrate lattice. But, as in the case of 
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Fig. 1. Titanium (001) surface with adatoms in a P(2x2) A-top 

configuration. Rectangles in this and subsequent figures 
indicate impact zones used in the simulation runs. 
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Fig. 2. Titanium with adatoms in a C(2x2) A-top configuration. 



29 



I 




Fig. 3. Titanium with adatoms in P(2x2) bridge configuration. 
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Fig. 4. Titanium with adatoms in C(2x2) bridge configuration. 
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Fig. 5, Titanium with adatoms in P(2x2) three-fold "A" position. 
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Fig. 6. Titanium with adatoms in C(2x2) three-fold "A" position. 
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Fig. 7 




. Titanium with adatoms in P(2x2) three-fold "B" position. 
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Fig. 8 




1 



. Titanium with adatoms in C(2x2) three-fold "B" position. Small 
dots indicate centers of atoms in the second layer. 



35 



Fig. 9 




Vanadium and Niobium clean (110) surface. 



36 




Fig. 10. Vanadium and niobium with adatoms in P(2x2) A-top position. 
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Fig. 11. Vanadium and niobium with adatoms in C(2x2) A-top position. 
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TABLE 2 



Physical Properties of Vanadium and Parameters 
of Associated Potential Functions 



Atomic Number 




23 


REF 


Atomic Weight 




50.9415 


30 


Melting Point 




1890 ± 10° C 


30 


Boiling Point 




3380° C 


30 


Valance 




+2, 3, 4, 5 


30 


Crystal Type 




bcc „ 


31 


a 




3.0399 ± 0.0003 A 




Lattice Unit 




1.51995 A 


Calc . 


^e 




2.6326 A 


31 


Cohesive Energy 




5.31 eV 


33 


Gas Phase Bond 


O 

Lengths (A) 






V-V 




2.6224 


30 


Ti-0 




1.56 ± .04 


34 


0-0 




1.20752 


35 


Ar-0 




3.31 


35 


Gas Phase Bond 


Strengths (eV) 






V-V 




2.511 


30 


V-0 




6.462 


30 


0-0 




5.1545 


35 


Ar-0 




0.0049 


30 


Surface Binding 


Energy 






V-0 




5.52 eV 


Calculated 



see text. 
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TABLE 3 



Physical Properties of Niobium and Parameters 
of Associated Potential Functions 



Atomic Number 


41 


REF 


Atomic Weight 


92.9064 


30 


Melting Point 


2468 ± 10® C 


30 


Boiling Point 


4742® C 


30 


Valance 


+2,3,4(?) ,5 


30 


Crystal Type 
a 


bcc o 

3.30656 ± 0.00002 A 


31 


Lattic Unit 


1.65328 A 


Calc 


^e 


2.8636 A 


Calc 


Cohesive Energy 


7.57 eV 


33 


O 

Gas Phase Bond Lengths (A) 


Nb-Nb 


2.8584 


30 


Nb-0 


1.69 


34 


0-0 


1.20752 


35 


Ar-0 


3.31 


35 


Gas Phase Bond Strengths (,eV) 


Nb-Nb 


2.83 


34 


Nb-0 


7.805 


35 


0-0 


5.1545 


35 


Ar-0 


0.0049 


30 



Surface Binding Energy 

Nb-0 6.67 eV Calculated, 

see text. 
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titanium, the reference index for the adatom mesh is uncer- 
tain. Two locations are possible on the vanadium and niobium 
(110) faces; the A- top and a rectangular four-fold site not 
present in titanium. In the four-fold position the adatoms 
could be placed either equidistant from the two pairs of 
neighbor atoms, or off to one side in a three-fold 'plus' 
position. Both would result in a high degree of bonding 
asymmetry. Since there is no corresponding site on the Ti 
(001) surface, the A- top position was chosen for comparison 
with titanium. 

No information is available on the binding energies of 
chemisorbed oxygen on vanadium or niobium. Since all three 
metals possess similar chemical properties, surface binding 
energies were calculated by scaling the gas phase M-0 
binding energies to the known SBE of oxygen on titanium. 

This resulted in estimated SBE's of 5.52 and 6.67 eV 
respectively. 

D . ADATOM COVERAGE 

Computations were run at two specific coverages, P(2x2) 
and C(2x2), for each of the adatom locations. These cover- 
ages correspond roughly to a quarter and a half monolayer 
of oxygen respectively. 

The designations for surface coverages were standardized 
in an article by E.A. Wood in 1964 [Ref. 45J. In this nota- 
tion the three dimensional unit cell is reduced to a diperiodic 
structure referred to as a net in which the unit area is the 
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unit mesh. The fourteen Bravais space lattices are thus 
reduced to five basic nets pictured below. 





Obi ique 

Y /s 90’, a b 



T 1 


) c 


5 f 




b 


o 


Y 




Y / 


6 


) c 


^ 



a 

Rectangl e 
Y = 90’, a b 



a 

Centered Rectangle 
Y = 90’, a / b 



The letter "P" is used to denote the primative mesh with 
equivalent points at the four corners, but none within the 
mesh itself. The letter "C is used to identify a unit mesh 
that contains a central point. Numerals are used to indi- 
cate the repetition interval in the X and Z directions, 
referenced to the substrate unit mesh. Meshes rotated from 
the substrate principal directions are designated as above 
followed by the angle of rotation. NaCl (100) PC5 x1)l- 45 
would represent a primative mesh located on the (100). 
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surface of NaCl, five units by one unit in area and oriented 
in the <110> direction. Meshes formed on other than sub- 
strate atoms are further identified by their origin point, 
two-fold bridge, etc. (For simplicity the four number 
designation normally used for hep crystals has been abbre- 
viated to three digits; i.e., (0001) = (001). The numbers 

themselves refer to the principal crystal directions (hkil) 
where i = h + k.) 

In the case of the titanium (001) surface and the (110) 
surfaces of vanadium and niobium, the unit mesh is a 120 
degree rhombus. The unit mesh for the A- top positions is 
centered above the substrate atoms and is displaced for the 
two and three-fold positions. Figures 1-11 depict the various 
coverages used, with the adatoms scaled according to their 
covalent radii. The two-fold locations are abbreviated "BR” 
(bridge) and the three-fold sites by ''tf”. Since there are 
two different three-fold sites on the hep (001) surface, the 
additional identifier "TFA" or "TFB” is used to denote 
whether the adatom is fixed "above” an atom in the second 
layer of the substrate or centered over a channel . 

E. POTENTIAL FUNCTIONS 

1 . General 

A variety of potential function forms are available 
in this program to describe the interatomic forces. The 
merits of each have been hotly debated in the literature, 
but Harrison has found [Ref. 42] that the global results are 
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not strongly dependent on the particular form in use. Indi- 
vidual parameters may be adjusted, if necessary, so that 
absolute values match experimental results, but the basic 
mechanisms involved remain unchanged. The four potential 
functions used are as follows: 

a. Born-Mayer , which is strictly repulsive having the 
form 



V(r) = A exp(-br) ; 

b. Moliere , which is an approximation to the Thomas-Fermi 
screening function, and is also repulsive 

V(r) = (Z Ze^/r) {0.35 exp (-0.3r/a) + 0 . 55 exp (-1 . 2r/a) 

1 ^ 

+ 0 . 1 exp (-6 . Or/a) } , 

where: 

a = 0.8853aQ/(z]^^^ + (Firsov radius) 



and 



= 0.5292 A. (Bohr radius) 

The parameter "a” may be varied to produce a "modified" 
Moliere function. 
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c. Mor se , which is both rspulsive and attractive having 
the form 



V(r) - D exp {-2a (r-r ) } - 2D exp{-a(r-r )} / 

where r^ is the equilibrivim separation distance between 
atoms, is the well depth and alpha is used to scale the 
distance in the potential function; and 

d. Composite Mol ier e-Morse , which consists of the two 
potential functions joined together smoothly by a cubic 
spline. 

The adjustable parameters can be related to the ion 
"size" and a collision "hardness", which in the case of the 
Born-Mayer potential correspond to the A and b constants. 

Thus one can tailor the frequency and hardness of collisions 
by appropriate selection of the parameters. At large separa- 
tions the Moliere and Born-Mayer functions have similar shapes. 
In all of the potentials the repulsive wall dominates the 
dynamics of collisions, while the attractive well determines 
whether an atom will escape the svirface and be sputtered. 
Reference 142] contains a detailed discussion on potential 
function selection. 

2 . Solid Phase Potentials 

The Ti-Ti interactions within the crystal were 

O 

modeled using a standard Moliere-Morse spline (a = 0.10529 A) 
to allow for the anticipated hard collisions in the experi- 

O 

ment. R^ was set at the lattice parameter value of 2.95 A 
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ad was adjusted to give the correct crystal cohesive 
eergy of 4.855 eV. The Morse alpha was set at 1.33 
S' that the slope closely matched that of the Moliere. The 
tfo were then connected by a spline between R = 0.68 and 
fj = 0.72 LU. the distance at which the potential is 

truncated/ was set at 2.20 LU to allow only interactions be— 
^^?een nearest and next-nearest neighbor atoms. Since the 
orces produced at greater separations are very weak, this 
election of R^ allows calculation time to be minimized with 
egligible effect on the results. 

The Nb-Nb and V-V potentials were similarly chosen 

0-1 0^1 

ath alpha equal to 1.34A and 1.62 A respectively. R^^ 

ras set at 2.2 LU in both cases. R^ was calculated as /3/2 

e 

:imes the lattice constant for the bcc structure. 

The metal (M) -Ar potential functions used for the 
Oombardment of the crystal by argon were chosen as standard 
Moliere potentials since no binding is to be expected. R^ 
was set at 1.7 LU for all runs. 

Morse potentials were used for the M-0 surface 
interactions. Lacking any specific data in the literature, 
adatoms in the A-top position were placed at the M-0 gas 
phase equilibrium separation distance above the substrate, 
although this is known to be too high. The well depth and 
alpha were then adjusted so that the correct surface binding 
energy was obtained. For the titanium two-fold bridge, gas 
phase equilibrium separation was maintained between the oxygen 
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and its two neighbor titaniums which necessitated placing the 
oxygens 0.322 LU above the surface plane. The well depth 
was again adjusted for the correct SBE. Adatoms in the 
three-fold sites were placed in the center of the titanium 
triangles on the surface plane, which required increasing 

O 

r to 1.70 A. R„ for all was set just outside the nearest 
neighbor distance to prevent 0-0 interference. Parameter 
values for the various adatom positions are summarized in 
Table 4 . 

The Ar-0 potential used is a Born-Mayer with A = 14.67 
o 

keV and b = -4.593 A . This choice gave good results in 

previous work by Harrison, et al. [Ref. 21J . 

The 0-0 potential is another Morse with = 5.1545 

® — 1 

eV, r^ = 1.208 A and alpha = 2.691 A ; values which were 
obtained from gas phase data and were also used by Harrison, 
et al . 

TABLE 4 

. Solid Phase Potential Function Parameters 







D (eV) 
e 




Alpha (A”^i 


Rc (1). 


Ti 


ATOP 


5.90 


1.62 


2.0 


2.086 


Ti 


BRIDGE 


2.95 


1.62 


2.0 


2.086 


Ti 


TF 


1.97 


1.70 


2.0 


2.086 


V 


ATOP 


5.52 


1.56 


2.0 


2.584 


Nb 


ATOP 


6.67 


1.69 


2.0 


2.810 
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3 . Gas Phase Potentials 

The potentials used for interactions above the 
surface, as well as in determining multimer formation, are 
of the pure Morse type. Parameters were selected from 
References 27, 32 and 33 and are summarized in Table 5. 
Appendix A (Figs. 15-22), contain plots of potential energy 
versus separation for most of the potential functions. 



TABLE 5 

Gas Phase Potential Function Parameters 



O 0—1 





Dg (eV) 


(A) 


Alpha (A 


Ti-Ti 


1.466 


2.896 


1.330 


Ti-0 


6.870 


1.620 


2.000 


Nb-Nb 


2.830 


2.858 


1.340 


Nb-0 


7.805 


1.690 


2.000 


V-V 


2.511 


2.622 


1.620 


V-0 


6.462 


1.560 


2.000 


0-0 


5.155 


1.208 


2.691 



F . MICROCRYSTAL 

Microcrystals composed of 676 atoms were used for each 
computer run. The atoms are placed by a lattice generator 
into the proper hep or bcc arrangement, resulting in a 
crystal 15 planes wide by 15 planes long by six planes deep 
(15x6x15) . The first atom in the crystal is number two, 
one being the ion, and is located at (0,0, 0,0). Atoms 2-114 
comprise the top layer, 115-226 layer two, 227-339 layer three. 
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340-451 layer four, 452-564 layer five and 565-676 layer 
six. Atom numbers above 676 correspond to oxygen adatoms. 

Crystal size selection is based on the problem of confine- 
ment which is the attempt to encompass the entire surface 
sensitive portion of the collision cascade within the target 
lattice. As the target size is increased the yield will 
asymptotically approach a maximum value. Past this point 
further increases in size will only affect computer run times. 
The 15x6x15 crystal size was selected after making a series 
of runs on clean titanium using a high yield impact point. 
Yield versus ion energy and crystal size is depicted in 
Figxares 12 and 13. As can be seen containment is effectively 
achieved for ion energies up to 3.0 keV. 

G. IMPACT AREA 

The impact area for each simulation is a rectangular 
array of 104 uniformly distributed impact points adapted from 
previous work [Ref. 46] . Figure 14 is a graphical represen- 
tation of the points and their coordinates. The actual length 
and height of the grid is automatically scaled by the pro- 
gram so that its dimensions match those of the interatomic 
spacing in the X and Z directions. Thus a grid positioned 
on planes RBX = 6.0 and RBZ = 6.0 would have its four corners 
located on the centers of the atoms at (6,6), (8,6), (6,8) 

and (8,8) . 

The impact points are slightly offset in the X and Z 
directions from true center to prevent repetitive sampling 
of symmetrically equivalent points. 
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Fig. 12, Yield versus crystal size at various ion energies. 
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Fig. 13 Yield versus crystal depth. 
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Fig. 14 



Standard impact zone with 104 



impact points. 
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Independent data sets can be generated by simply moving 
the origin of the grid. Selection of grid position is based 
on the surface symmetry resulting from adatom density and 
placement. Generally two sets of impact points are re- 
quired to sample a completely representative area of the 
target surface. Surfaces which could be mapped into one 
grid were still run with two sets, shifted by one plane, 
so that all runs consisted of 208 impact points. 

H. ION ENERGY 

All simulations were conducted with an ion energy equal 
to 1.0 keV. This energy produced a reasonable sputtering 
yield and allowed the runs to be conducted on a moderate 
sized crystal. Additional runs were conducted on the 
titanium C(2x2)TFB surface at 0.5 and 2.0 keV, for compari- 
son. This coverage produced the highest yield at 1.0 keV 
and therefore the differences in results as the ion energy 
is varied should be more pronounced. 

I. PROGRAM MODIFICATIONS 

1. QRAD 

During the early stages of this investigation 
Harrison completed extensive modifications to QRAD resulting 
in improved energy conservation, timestep logic and output 
format. Considerable effort was expended proofing and de- - 
bugging the new version. Additional code was also inserted 
to identify registered pairs, metal-oxygen pairs on the 
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undamaged surface, for later use in tracking multimer 
origins. 

2. AN83 

This is the analyzer program used to determine vari- 
ous yield distributions and multimer data. Changes were 
made to make the program compatible with the new output 
format of QRAD. 

3. PLOT AN 

This program used the output from QRAD to generate 
graphic displays of the results. The program was modified 
to permit separate analyses of each atom type. It was also 
expanded to accommodate additional crystal types. The logic 
was modified to allow determination of the angular distributions 
from data generated by the non-symmetric sets of impact 
points in use. Plot scaling routines were also changed to 
handle the additional crystal types. 

4 . ADPLOT 

ADPLOT is used to generate the pictorial representa- 
tion of adatoms and steps on the crystal surface. Modifica- 
tions were made to permit real time changes to the scaling 
of the plot size and to adjust adatom radii. 
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IV. RESULTS 



A. GENERAL 

Detailed information on each run has been assembled into 
appendices following the main body of this thesis. Included 
are graphical displays of the ejection time distributions, 
atoms ejected per single ion (ASI) , energy distributions, 
spot patterns, and the ejection time differences between 
multimer components. Where applicable these have been broken 
down into separate substrate and adatom plots. Plot headings 
are generally self-explanatory, indicating ion energy, metal 
type and crystal face used for the particular run. For 
example "1.0 keV Ti(001)/Ar + Ox P(2x2)" indicates that the 
data pertains to bombardment of the 001 face of titanium 
by 1.0 keV argon ions. Ox PC2x2) denotes the presence of 
oxygen adatoms in a PC2x2) configuration. (M) or (OX) are 
used to identify whether the plot is for the substrate or 
for the adatoms in the particular configuration. The plots 
contain information based on a standard run of 208 impact 
points. On occasion data from a multiple set of 520 impact 
points is included to amplify significant characteristics. 
These are identified by the notation (IPX51 . Combined plots 
of adatoms in the TFA and TFB positions are identified by 
TF CA+B) . Discussions of the more noteworthy results of the 
experiments are contained in the remainder of this chapter. 
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B. YIELD AND AS I 



The yields from the runs conducted on the clean crystal 
surfaces are summarized in Tables 6 and 7. All runs consisted 
of the standard set of 208 impact points with the exception 
of the expanded run on titanium which was composed of 520 
impact points. This larger run was used to test the sta- 
tistical validity of findings from the smaller sets. Of 
particular interest is the fact that the yield from the 
representative sample of 208 impact points for titanium 
is within 1.6% of the yield from a much larger run of 520 
impact points, demonstrating the validity of this rather small 
set. The significantly higher yield of argon from the bom- 
bardment of niobium and vanadixm is to be expected in view 
of the lower mass ratio, 0.429 Cversus 0.834 and 0.784 for 
Ti and V respectively) and the increase in the collision 
cross-section which is a function of the atomic number. 

The increased yield of V and Nb, as compared to Ti, may 
be explained by two characteristics; 

1. The more open nature of the bcc (.110) face. While 
the overall cohesive energy is greater, the surface layer 
has fewer nearest neighbors and therefore fewer bonds must 
be broken for an atom to escape the surface. 

2. The fact that there are no channels present in the 

bcc CllO) face for bombardment at normal angles of incidence. 
Hence all the ion energy is deposited in the top layers of 
the crystal from which the vast majority of sputtered atoms 
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TABLE 6 



Yield and ASI (by Coverage) 



COVERAGE 

e 




(%) 


^0 


(%) 


Nar 




^0 


CLEAN 
















Ti (520) 


945 








10 


1.8 




Ti (208) 


374 




— 




2 


1.8 


_ 


V 


946 




- 




8 


4.5 


— 


Nb 


682 




- 




52 


3.3 


- 


P(2x2) 
















Ti ATOP 


291 


(79) 


77 


(21) 


2 


1.4 


0.4 


V ATOP 


616 


(83) 


126 


(17) 


6 


3.0 


0.6 


Nb ATOP 


522 


(82) 


112 


(18) 


37 


2.5 


0.5 


Ti BR 


331 


(80) 


82 


(20) 


4 


1.6 


0.4 


Ti TEA 


324 


(90) 


38 


(10) 


2 


1.6 


0.2 


Ti TFB 


386 


(87) 


58 


(13) 


4 


1.9 


0.3 


C(2x2) 
















Ti ATOP 


215 


(62) 


133 


(38) 


2 


1.0 


0.6 


V ATOP 


457 


(70) 


198 


(30) 


6 


2.2 


0.9 


Nb ATOP 


374 


(70) 


159 


(30) 


41 


1.8 


0.7 


Ti Br 


271 


(65) 


14 8 


(35) 


4 


1.3 


0.7 


Ti TEA 


290 


(77) 


86 


(23) 


4 


1.4 


0.4 


Ti TEB(0.5) 


233 


(69) 


105 


(31) 


7 


1.1 


0.5 


Ti TEB(l.O) 


368 


(75) 


124 


(25) 


6 


1.8 


0.6 


Ti TEB(2.0) 


419 


(78) 


116 


(22) 


2 


2.1 


0.6 



NOTE: All ion energies are 1.0 keV unless otherwise noted, 

i.e., (0.5). The designations (208) and (520) 

refer to a single set of impact points and a 
combined group of 520 impact points respectively. 

N, is the number of reflected argon ions. N„ 

and Nq are the actual counts of ejected atoms. 

and Yq are the computed average yields, atoms/ ion. 
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TABLE 7 



Yield and ASI (by Sxibstrate) 



COVERAGE 


Nm 


(%) 


1 


[%) 


N 


Y 


Y 




M 




0 




AR 


M 


0 


CLEAN 
















Ti (520) 


945 








10 


1.8 





Ti (208) 


374 




- 




2 


1.8 


- 


P (2X2) 
















Ti ATOP 


291 


(79) 


77 


(21) 


2 


1.4 


0.4 


Ti BR 


331 


(80) 


82 


(20) 


4 


1.6 


0.4 


Ti TF <A+B> 


355 


(88) 


48 


(12) 


3 


1.8 


0.3 


CC2X2) 
















Ti ATOP 


215 


(62) 


133 


(38) 


2 


1.0 


0.6 


Ti BR 


271 


(65) 


148 


(35) 


4 


1.3 


0.7 


Ti TF <A+B> 


329 


(76) 


105 


(24) 


5 


1.6 


0.5 


V CLEAN 


946 









8 


4.5 


_ 


V P(2X2) ATOP 


616 


(83) 


126 


(17) 


6 


3.0 


0.6 


V C(2X2) ATOP 


457 


(70) 


198 


(30) 


6 


2.2 


0.9 


Nb CLEAN 


682 




— 




52 


3.3 


— 


Nb P(2X2) ATOP 


522 


(82) 


112 


(18) 


37 


2.5 


0.5 


Nb C(2X2) ATOP 


374 


(70) 


159 


C30) 


41 


1.8 


0.7 



NOTE: The quantity <A+B> for the three-fold sites 

refers to the coinbined average of the A and B 
positions of the adatoms. 
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The yield per layer data for adatoms in the A- 
top position is given in Table 8. Yield ratios for the 
other adatoms locations were found to be essentially identical. 



TABLE 8 



Sputtering Yield by Layer 



COVERAGE 


LAYER 1 


LAYER 2 


LAYER 3 


LAYER 4 


CLEAN 










Ti (001) 


372 


0 


0 


0 


V (110) 


900 


44 


8 


0 


NB (110) 


694 


26 


8 


0 


P(2x2) 










Ti (001) 


285 


5 


1 


0 


V (110) 


566 


41 


7 


2 


Nb (110) 


490 


24 


8 


0 


C (2x2) 










Ti (001) 


213 


2 


0 


0 


V (110) . 


429 


24 


8 


1 


Nb (110) 


34 7 


24 


3 


0 


With respect 


to the various adatom 


coverages the 


metal 


yields for the A 


-top positions are fairly consistent 


at 


about 80% for the P(2x2) and 


70% for the C(2x2). As 


the 


oxygen coverage 


is increased 


there are 


fewer direct 


impacts 



on the substrate atoms which begin to look more like a second 
layer to the incoming ions with a corresponding decrease in 
yield. This decrease is summarized below where Delta Yj^ 
represents the drop in substrate yield as the coverage is 
increased from a quarter to a half monolayer. 
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TABLE 9 



Change in Metal Yield as a Function of Adatom Coverage 



COVERAGE 



DELTA (%) 



Ti P C ATOP 
V V C ATOP 
Nb P C ATOP 



-17 

-13 

-12 



Ti P C BR 

Ti P ^ C TFA 
Ti P C TFB 



-15 

-13 

-12 



As the adatom position is shifted from A-top to bridge 
to three-fold, the total yield (M+0) fluctuates with the 
maximum coming from the bridge positions. The metal yields 
are about the same percentage for both the A-top and bridge, 
but increase in the three-fold positions. The corresponding 
decrease in the fraction of oxygens ejected is reasonable 
in light of the added stability of the triply bonded positions. 

The increased metal yield in the bridge, and particularly 
in the TFB positions, can again be explained in terms of 
channel blockage. The bridge locations partially obstruct 
the channel in the hep CO 01) plane, while the TFB position 
above the channel closes it to incoming ions. 

The Atoms per Single Ion graphs for selected cases are 
contained in Appendix B. The distributions for the substrate 
atoms are in general broader than for the adatoms. As the 
oxygen density is increased, there is little change in the 
substrate ASI distribution, but the maximum in the oxygen 
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distributions shifts toward an ASI of two. This is to be 
expected due to the larger number of oxygen atoms on the 
surface. The ASI distributions of the substrate broaden 
somewhat as the adatom positions are changed from the bridge 
to the three-fold, consistent with the increasing metal 
yield. 

C. EJECTED PARTICLE ENERGY VERSUS COVERAGE 

It is interesting to note that as the adatom coverage 
in the A-top position is increased, the average energy of 
the ejected substrate atoms also increases. As noted above, 
this corresponds to a decrease in metal atom yield as well. 
Lower energy particles are apparently retained in the crystal, 
thus raising the average ejection energy. The most probable 
energy of the ejected metal atoms also is significantly 
higher when there are adatoms in the A-top position and in- 
creases with oxygen density as well (Table 10) . This trend 
does not extend to the bridge and three-fold positions in 
titaniiam, but the most probable energies of the substrate 
atoms are more in line with those ejected from clean metal 
surfaces . 

There is a much closer correlation between the metal and 
oxygen average energies in the bridge and three-fold posi- 
tions. The multiple bonding involved appears to equalize 
the energies of the ejected particles to some degree. 

As the ion energy is increased there is a corresponding 
rise in the average ejection energy of the metal atoms, as 
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TABLE 10 



Average and Most Probable Ejection Energies 



COVERAGE 


<Ej^> (eV) 


[E^] MP (eV) 


A 

W 

o 

V 


CLEAN 








Ti 


10.3 


3.5 




V 


8.9 


4.0 


— 


Nb 


12.3 


4.5 


- 


P(2x2) (ATOP) 








Ti 


10.7 


7.0 


11.9 


V 


10.7 


7.5 


11.4 


Nb 


13.8 


9.0 


11.9 


C(2x2) (ATOP) 








Ti 


11.3 


7.5 


27.5 


V 


12.2 


8.0 


10.7 


Nb 


15.9 


9.0 


13.8 


TITANIUM 








P(2x2) BR 


10.057 


4.0 


10.062 


C(2x2) BR 


9.815 


5.0 


9.461 


P(2x2) TEA 


9.852 


5.0 


10.997 


C(2x2) TEA 


10.193 


4.5 


10.190 


P(2x2) TFB 


10.230 


3.0 


8.055 


C(2x2) TEB 


9.743 


3.0 


9 .384 
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TABLE 10 



Average and Most Probable Ejection Energies 



COVERAGE 


<Ej^> (eV) 


[E^] MP (eV) 


A 

W 

o 

V 


CLEAN 








Ti 


10.3 


3.5 





V 


8.9 


4.0 


— 


Nb 


12.3 


4.5 


- 


P(2x2) (ATOP) 








Ti 


10.7 


7.0 


11.9 


V 


10.7 


7.5 


11.4 


Nb 


13.8 


9.0 


11.9 


C(2x2:) (ATOP) 








Ti 


11.3 


7.5 


27.5 


V 


12.2 


8.0 


10.7 


Nb 


15.9 


9.0 


13.8 


TITANIUM 








P(2x2) BR 


10.057 


4.0 


10 .062 


C(2x2) BR 


9.815 


5.0 


9.461 


P(2x2) TEA 


9.852 


5.0 


10.997 


C(2x2) TEA 


10.193 


4.5 


10.190 


P(2x2) TEB 


10.230 


3.0 


8.055 


C(2x2) TEB 


9.74 3 


3.0 


9.384 



62 



TABLE 11 



Average Ejection Energy Versus Ion Energy 



ION ENERGY (keV) 



Ti (eV) 



Ox (eV) 



0.5 

1.0 

2.0 



8.0 59 
9.743 
13.146 



22.514 

9.348 

20.316 



one would expect. The changes in the oxygen energy match 
the variation in yield (105, 124, 116) . While the ratio of 
energies between Ti and 0 is essentially the same for 0.5 
and 2.0 keV, the oxygen energy is about double that for the 
1.0 keV ions at all coverages except the A-top C(.2x2) where 
there is also an unexpected rise. 

The energy distribution curves are contained in Appendix 

C. The curves for the ejected oxygen atoms are typically 
flat for the bridge and the three-fold positions, showing 
no dominant ejection energy. However there is a pronounced 
Rayleigh-type shape to the curves for oxygen in the A-top 
position and they bear more of a resemblance to those of the 
substrate. This agrees with what was found earlier where 
the multiple bonds tend to cause a more even distribution of 
ejected particle energy. 

D. EJECTION TIME DISTRIBUTIONS 

Appendix D contains bar charts of the ejection time dis- 
tributions for representative runs . The distributions are 
similar for titanium and vanadixam, but the average ejection 
time for niobium increases about 30 fsec to about 110 fsec 
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sec) . The distribution is also 



-15 

(1 femptosecond = 10 
broader. There is essentially no change in any of the dis- 
tributions as oxygen is added to the surface, or increased 
from P(2x2) to C{2x2). Thus the ejection time appears to 
be mainly a function of the substrate mass and independent 
of adatom coverage and lattice type. 

The oxygen ejection time distributions show the same 
characteristics as the substrate atoms. The only significant 
difference was that oxygen ejection started approximately 
10 fsec later for titanium than for vanadium or niobixom. 

Again the distribution was broader for niobium, reflecting 
the influence of the substrate. 

E. YIELD VERSUS ENERGY 

Since the C{2x2) TFB coverage produced the highest titanium 
sputtering yield, it was chosen to investigate the effects 
of changing ion energy. This coverage was therefore run 
at energies of 0.5 and 2.0 keV in addition to the 1.0 keV 
trial . The yield results and layer data are contained in 
Tables 12 and 13. 



ENERGY 

(keV) 

0.5 

1.0 

2.0 



TABLE 12 

Titanium Atoms Ejected and Yield 
as a Function of Ion Energy 




M 



N, 



233 

368 

419 



105 0.5 
124 0.6 
116 0.6 



N 

0 

0 

0 



Ar 



1.1 

1.8 

2.0 
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TABLE 13 



Titanium Atom Yield per Layer 
as a Function of Ion Energy 



ENERGY 

(keV) 


LAYER 1 


LAYER 2 


LAYER 3 


LAYER 4 


0.5 


232 


1 


0 


0 


1.0 


363 


4 


1 


0 


ro 

• 

o 


39 8 


16 


1 


4 



As can be seen the higher the energy, the higher the yield, 
but the relationship is definitely non-linear. With increasing 
energy the average ion penetrates further into the crystal . 

As a consequence a smaller percentage of its energy is deposited 
in the top layers. More energy is carried deep into the 
crystal where it is dissipated through collision cascades 
producing only lattice displacements and vibrations. 

The decrease in oxygen yield for the 2.0 keV ions is 
also probably caused by the deeper deposition of energy. 

The energy distribution plot for the ejected oxygens (Figs. 
54-56) clearly shows a decrease in the energy as compared to 
the 1.0 keV case. Hence fewer adatoms are receiving enough 
energy to surmount the surface binding well and escape. 

F. MULTIMER YIELD 

Tables 14 and 15 present the various multimer yields in 
terms of total numbers and as percentages. Tables 16 and 17 
contain a breakdown by species arranged by substrate type 
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TABLE 14 



Multimer Yield as a Function of Coverage 



COVERAGE 


1-MER 


2-MER 


3-MER 


4 -HER 


5-MER 


TOTAL 

MULTIMER 


Clean 


Ti (520) 


860 


43 


2 


0 


0 


45- 


Ti (208) 


356 


10 


2 


0 


0 


12 


V 


731 


74 


22 


2 


0 


98 


Nb 


585 


48 


13 


2 


0 


63 


P(2x2) 


Ti ATOP 


304 


25 


4 


0 


0 


29 


V ATOP 


529 


77 


12 


5 


0 


94 


Nb ATOP 


441 


60 


18 


2 


0 


80 


Ti BR 


350 


26 


3 


0 


0 


29 


Ti TFA 


337 


10 


1 


0 


0 


11 


Ti TFB 


389 


23 


1 


2 


0 


26 


C(2x2) 


Ti ATOP 


250 


49 


2 


0 


0 


51 


V ATOP 


416 


99 


3 


3 


2 


107 


Nb ATOP 


377 


59 


7 


2 


0 


68 


Ti BR 


335 


31 


3 


2 


0 


36 


Ti TFA 


336 


16 


2 


0 


0 


18 


Ti TFB ( .5) 


310 


14 


0 


0 


0 


14 


Ti TFB (1.) 


437 


18 


5 


0 


0 


23 


Ti TFB (2.) 


490 


17 


3 


0 


0 


20 



NOTE: All ion energies are 1.0 keV except as indicated, 

i .e . , (0.5). 
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TABLE 15 



Multimer Yield (Percentage) 



% YIELD 
THAT FORMED 



COVERAGE 


1-MER 


2-MER 


3-MER 


4-MER 


5-MER 


MULTIMERS 


Clean 


Ti (520) 


90.4 


9.0 


0.6 


0.0 


0.0 


9.6 


Ti (208) 


97.9 


5.3 


1.6 


0.0 


0.0 


6.9 


V 


76.9 


15.4 


6.9 


0.8 


0.0 


23.1 


Nb 


80.5 


13.1 


5.3 


1.1 


0.0 


19.5 


P(2x2) 


Ti ATOP 


83.2 


13.5 


3.3 


0.0 


0.0 


16.8 


V ATOP 


71.6 


20.8 


4.9 


2.7 


0.0 


28.4 


Nb ATOP 


71.3 


18.9 


8.5 


1.3 


0.0 


28.7 


Ti BR 


85.2 


12.6 


2.2 


0.0 


0.0 


14.8 


Ti TFA 


93.7 


5.5 


0.8 


0.0 


0.0 


6.3 


Ti TFB 


87.1 


10.4 


0.7 


1.8 


0.0 


12.9 


C(2x2) 


Ti ATOP 


70.1 


28.2 


1.7 


0.0 


0.0 


29.9 


V ATOP 


65.1 


30 .2 


1.4 


1.8 


1.5 


34.9 


Nb ATOP 


72.5 


22.1 


3.9 


1.5 


0.0 


27.5 


Ti BR 


82.2 


14.8 


2.1 


1.9 


0.0 


18.8 


Ti TFA 


89 .9 


8.5 


1.6 


0.0 


0.0 


10.1 


Ti TFB ( .5) 


91.7 


8.3 


0.0 


0.0 


0.0 


8.3 


Ti TRB (1.) 


89.6 


7.3 


3.1 


0.0 


0.0 


10.4 


Ti TFB (2.) 


93. 6 


6.4 


1.7 


0.0 


0.0 


8.1 



NOTE: All ion energies are 1.0 keV except as indicated, 

i .e . , (0.5). 
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TABLE 16 



Multimer Species Formed (by Si;ibstrate) 



COVERAGE 



M, 




MO M 2 O M^O MO 2 



M2O2 






Clean 



Ti (520) 


43 


2 


0 






_ 






Ti (208) 


10 


2 


0 


- 


- 


- 


- 


- 


(2x2) 


Ti ATOP 


3 


0 


0 


22 


4 


0 


0 


0 


Ti BR 


8 


0 


0 


18 


3 


0 


0 


0 


Ti TF <A+B> 


7 


0 


0 


8 


1 


1 


0 


0 


(2x2) 


Ti ATOP 


0 


0 


0 


49 


2 


0 


0 


0 


Ti BR 


5 


0 


0 


26 


2 


1 


1 


1 


Ti TF <A+B> 


6 


1 


0 


11 


2 


0 


1 


0 


V CLEAN 


74 


22 


2 













V P(2x2) 


18 


0 


2 


59 


12 


3 


0 


0 


V C(2x2) 


9 


0 


0 


90 


3 


1 


0 


2 


Nb CLEAN 


48 


13 


2 











_ 


Nb P(2x2) 


18 


6 


0 


42 


12 


1 


0 


1 


NB C(2x2) 


2 


2 


0 


57 


5 


0 


1 


2 



0 

2 



0 

0 



NOTE: The quantity <A+B> refers to the average of the two 

three- fold positions. 
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TABLE 17 



COVERAGE 



Clean 

Ti (520) 
Ti (208) 

V 

Nb 

P (2x2) 

Ti ATOP 

V ATOP 
Nb ATOP 
Ti BR 
Ti BFA 
Ti TFB 

C(2x2) 

Ti ATOP 

V ATOP 
Nb ATOP 
Ti BR 
Ti TFA 
Ti TFB ( 
Ti TFB ( 
Ti TFB ( 



NOTE; 



Multimer Species Formed (by Coverage) 






M. 






iMO 



M^O 



M3O 



MO, 



M2O2 



^ 3°2 



43 2 0 

10 2 0 

74 22 2 

48 13 2 



3 0 

18 0 

18 6 

8 0 

6 0 

11 0 



0 0 

9 0 

2 2 

5 0 

6 0 

.5) 4 0 

L.) 7 1 

2.) 9 0 



0 22 4 

2 59 12 

0 42 12 

0 18 3 

0 4 1 

0 12 1 



0 49 2 

0 90 3 

0 57 5 

0 26 2 

0 10 2 

0 10 0 

0 11 2 

0 8 2 



0 0 

3 0 

1 0 

0 0 

0 0 

2 0 



0 0 
1 0 
0 1 
1 1 
0 0 
0 0 
0 2 
0 1 



All ion energies are 1.0 keV except as indicated, 
i.e. , (0.5) . 
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OOiHOOO OCNCNiHOOOO 



s.nd ddatom covsrags. The inuXtimer yield is seen to increase 
with oxygen coverage and is due primarily to the ejection 
of metal atoms with their attached adatoms. There is a good 
co^^®spondence in the mul timer yield between vanadium and 
niobium in the clean and P (2x2) coverage which then degener- 
ates as the coverage is increased to C(2x2) . This can be 
attributed to the decreasing influence of the crystal lattice 
and increasing effect. of the oxygen on the ejection mechan- 
isms. This will be covered in greater detail in the 
following section. 

G. MULTIMER ORIGINS 

One of the primary objectives was to study the origins 
of multimer components to compare fragmentation and recom- 
bination as modes of molecule formation. Table 18 presents 
data as to whether the constituent atoms were nearest-neighors 
in the crystal prior to ejection, or in the case of trimers 
or greater, whether a portion of the mdlecule was joined. 

In the table the letter 'C represents the number of multimers 
whose components were coupled together on the surface; p — 
partially connected multimers, N — those multimers whose 
components were not joined before ejection, %1 — the percentage 
that was formed by fragmentation and %2 — the percentage 
that were formed either through total or partial fragmentation, 
where a portion of the multimer was joined on the surface. 

Table 19 contains similar data for trimers and higher order 
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TABLE 18 



Fragmentation versus Recombination 
(All multimers) 



COVERAGE 


C 


P 


N 


%1 


%2 


Clean 


Ti (520) 


29 


1 


15 


64.4 


66.6 


Ti (208) 


6 


1 


5 


50.0 


58.3 


V 


71 


7 


20 


72.4 


79.6 


Nb 


41 


1 


21 


65.1 


66.7 


P(2x2) 


Ti ATOP 


24 


0 


5 


82.8 


82.8 


V ATOP 


76 


6 


12 


80.9 


87.2 


Nb ATOP 


52 


7 


21 


65.0 


73.8 


Ti BR 


23 


2 


4 


79.3 


86.2 


Ti TFA 


9 


1 


1 


81.8 


90.9 


Ti TFB 


21 


0 


5 


80.8 


80.8 


C(2x2) 


Ti ATOP 


49 


0 


2 


96.0 


96.0 


V ATOP 


86 


2 


19 


80 .4 


82.2 


Nb ATOP 


49 


4 


15 


72.1 


77.9 


Ti BR 


29 


2 


5 


80 .6 


86.1 


Ti TFA 


11 


0 


7 


61.6 


61.6 


Ti TRB (.5) 


10 


0 


4 


71.4 


71.4 


Ti TFB (1.) 


17 


1 


5 


73.9 


78.3 


Ti TFB (2.) 


13 


1 


6 


65.0 


70.0 



NOTE; C represents the number of multimers that were nearest 
neighbors on the surface prior to ejection; P, those 
partially connected and N, the mamber formed by 
recombination of non-adjacent atoms. %1 is the 
percent formed by fragmentation, %2 the percent formed 
by partial fragmentation. All ion energies are 1.0 
keV except as indicated, i.e., (0.5). 
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TABLE 19 



Fragmentation versus Recombination 
(Trimers and Greater) 



COVERAGE 


C 


P 


N 


%1 


%2 


Clean 


Ti (520) 


1 


1 


0 


50.0 


100.0 


Ti (208) 


1 


1 


0 


50.0 


100.0 


V 


17 


7 


0 


70 .8 


100.0 


Nb 


14 


1 


0 


93.3 


100.0 


P(2x2) 


Ti ATOP 


23 


0 


3 


88.5 


88.5 


V ATOP 


11 


6 


0 


64 .7 


100.0 


Nb ATOP 


13 


7 


0 


65.0 


100 .0 


Ti BR 


1 


2 


0 


33.3 


100.0 


Ti TFA 


0 


1 


0 


0.0 


100.0 


Ti TFB 


3 


0 


0 


100.0 


100.0 


C(2x2) 


Ti ATOP 


1 


0 


1 


50 .0 


50.0 


V ATOP 


6 


2 


0 


75.8 


100.0 


Nb ATOP 


5 


4 


0 


55.6 


100.0 


Ti BR 


4 


2 


0 


60.0 


100.0 


Ti TRA 


1 


0 


1 


50 .0 


50.0 


Ti TRB (1.) 


3 


1 


1 


60.0 


80.0 


Ti TFB (2.) 


1 


1 


1 


33.3 


66.7 


NOTE; C represents 


the number 


of 


multimers that were 


nearest 


neighbors on the 


! surface prior 


to ejec 



tion; P, those partially connected and N, the 
number formed by recombination of non-adjacent 
atoms . %1 is the percent formed by fragmenta- 

tion, %2 the percent formed by partial fragmenta 
tion. All ion energies are 1.0 keV except as 
indicated, i.e., (0.5). 
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niultiiners. Table 20 presents the percentage of multimer 
species formed by fragmentation. 

The general trend is toward more fragmentation as the 
oxygen coverage is increased. This is due primarily to the 
fact that substrate atoms are ejected after being hit from 
underneath and carry off their attached adatoms as they leave. 
Many partial fragmentation molecules were formed when attached 
atoms exited the surface obliquely, knocking off oxygens from 
nieghboring atoms on the way out. This also appears to be 
one of the principal mechanisms involved in recombination 
of non-a joining atoms. 

Figures 74-87 of Appendix E and Table 21 give the loca- 
tions and frequency of occurrence for all of the observed 
mechanisms involved in multimer formation. (In the table 
double entries refer to the a and b sections of the particu- 
lar drawing.) By far the greatest source of multimers is 
the dimer composed of a substrate atom and its attached 
adatom which is a fragmentation mechanism. Higher order 
multimers are a fairly rare occurrence, but over 65% are 
formed by fragmentation and virtually all involve at least 
partial fragmentation (Table 20) . 

As the oxygen coverage is increased, the monoxide is 
still the prevalent species formed, but there is a shift in 
the mechanisms that make up the rest of the yield and new 
mechanisms are added. In the A- top and bridge positions 
the percentage of multimers resulting from fragmentation 
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TABLE 20 



Percentage of Multimer Species 
Formed by Fragmentation 



COVERAGE 


M 

^^2 


M 3 


M 4 


MO 


M^O 


M 30 


M 02 


M 2 O 2 


M 302 


Clean 

Ti 


67 


50 
















V 


74 


73 


100 


— 


— 






_ 




Nb 


56 


92 


100 


- 


- 


- 


- 


- 


- 


P(2x2) 


Ti ATOP 


33 


— 


— 


82 


100 






_ 




V ATOP 


100 


- 


100 


80 


50 


100 


— 


— 




Nb ATOP 


83 


100 


- 


57 


50 


100 


— 


IP 


— 


Ti BR 


63 


- 


- 


89 


33 


— 


— 


— 


— 


Ti TFA 


100 


- 


- 


75 


IP 


— 


— 


— 


— 


Ti TFB 


64 


- 


- 


92 


100 


100 


— 


— 


— 


Ti TF(A+B) 


82 


- 


- 


84 


50 


100 


- 


- 


- 


C(2x2) 


Ti ATOP 


— 


— 


_ 


98 


50 


— 




— 


— 


V ATOP 


56 


- 


- 


81 


67 


IP 


- 


100 


100 


Nb ATOP 


100 


100 


- 


74 


40 


— 


X 


X 


- 


Ti BR 


40 


- 


- 


92 


100 


100 


Ip 


IP 


- 


Ti TRA 


33 


- 


- 


80 


50 


— 


- 


— 


- 


Ti TFB 


86 


X 


- 


73 


100 


— 


50 


- 


- 


Ti TR(A+B) 


60 


X 


- 


77 


75 


— 


50 


- 


- 



NOTE: A tack (-) indicates that the species was not formed, 

while "X" indicates that none were formed by frag- 
mentation. IP indicates that only one of the particu- 
> lar species was formed and that was through partial 

fragmentation. 
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TABLE 21 



Frequency of Occurrence of the Various 
Ejection Mechanisms 



COVERAGE 


1 


2 


3 


4 


dean 


Ti (520) 


20/8 


1 


2 


9/3 


Ti (208) 


3/2 


2 


2 


1 


V 


37 


13 


4 


9 


Nb 


19 


8 


10 


4 


P(2x2) 


Ti ATOP 


1 


1 


19 


1/2 


V ATOP 


6/0 


6 


1 




NB ATOP 


5/1 


3 


3 


3 


Ti BR 


4 


1 


1 


1 


Ti TEA 


1 


2 


1 


1 


Ti TFB 


1 


5 


1 


1 



5 6 7 8 9 10 11 12 



11 - _____ 

5 7 - 4224- 

6 2 8 - 2--2 



2 47/10 - - 4 2 3 

4 - 24/7 7 1 5 1 2 

112 2 11 --- 

3 1 - _____ 

11 11111 -- 



C(2x2) 

Ti ATOP 
V ATOP 
Nb ATOP 
Ti BR 
Ti TEA 
Ti TFB (.5) 
Ti TFB (1.) 
Ti TTB (2.) 



1 45 1 1 

3 19 1 

2 - - - 

1111 
1.3 1 1 

- 1/1 - 
1 1/1 1/0 2 
1 - 3/2 - 



2-74 

- 2 42 

24 2 1 

1/1 6/2 1 

1 5/2 

1 4/4 1/1 

- 5/1 - 



5 7 3 -- 

10 2 - 1 2 

11111 

11111 
1 - - - 1 



NOTE; The column headings refer to surface locations shown in 

Appendix E. Columns with double entries refer to the a and 
b sections of the individual figures. Ion energies are 1.0 
keV except as indicated, i.e., (.5) . 
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TABLE 21 (CONT.) 
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Ti TFB (2.) 
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COVERAGE 25 26 27 28 

P(2x2) 

Nb ATOP 12 11 
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increases, while it decreases for the three— fold positions. 
Those oxygens, sharing three bonds, tend to be left behind 
unless at least two of their nearest neighbors are also 
ejected. The greater oxygen density also presents a larger 
target to obliquely ejected fragments, increasing the 
probability of collision and subsequent capture. 

H. MULTIMER EJECTION TIMES 

Appendix F contains bar charts depicting the differences 

in ejection times between multimer components, both for 

metal-metal and metal-oxygen pairs. Essentially none of 

the pairs in any of the runs were ejected simultaneously, 

but all but about seven pairs were ejected within 100 fsec 

of each other and most within 50 fsec. This difference is 

about the same order of magnitude as the vibrational period 

-13 

of an equivalent diatomic molecule [Ref. 47J , 10 sec 
(100 fsec) , and is not necessarily an indication that the 
atoms were ever disjoined at some point in time. 

Comparison of the charts for the different oxygen cover- 
ages reveals a correspondence between ejection time differ- 
ence and the percentage of multimers that were formed by 
fragmentation. As one would expect there is an increasing 
time difference between component ejection as the niomber of 
molecules formed by recombination increases. 

I . SPOT PATTERNS 

The spot patterns produced by ion bombardment of the 
crystal surfaces are contained in Appendix G. The ejection 
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p^ttsirns coirr'sspond to ths dixsctions of ths doss— packsd. 
rows in the crystal, which for hep are oriented at 30, 90 
and 150 degrees from the <00 1> direction. For bcc the direc- 
tions are 0 and 90 degrees. These lead to a six- fold symmetry 
pattern in hep and a four-fold rectangular pattern in bcc. 

The patterns for vanadium and niobium are in quite good 
agreement with the expected results, but there is consider- 
able deviation in the case of titanium atoms. The spikes at 
90, 150 and 330 degrees are present, but there are extra 
spikes at 0, 60, 180 and 240 degrees, while those at 30 and 
210 degrees are missing. A variety of different runs were 
conducted with additional impact point origins and permutated 
impact point sets, but all resulted in the same pattern. At 
present there is no apparent reason for this discrepancy. 
Second layer contributions would tend to broaden the spikes 
somewhat, but not change the overall appearance. 

The patterns for the metals with adatoms in the A-top 
position are essentially the same as for those from clean 
surfaces, even though they have considerably more energy as 
was pointed out earlier. There is some blurring, presumably 
due to adatom interference, but it is very minor. This 
blurring becomes more pronounced as the oxygen density is 
increased from P to C. The features of the metal spot 
patterns begin to fade when adatoms are placed in the bridge 
positions and are completely absent for oxygen in the three- 
fold sites, leaving a uniform angular distribution. 
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Thsre is no marked tendency for oxygen to eject perpendicu— 
Isr to the surface in any of the locations, as was found by 
Hanson [Ref. 37] . The ejection angles were observed to in- 
crease slightly with increasing oxygen density, probably 
caused by the interference of other adatoms and increasing 
yield. This was also found to be the case for varying ion 
energy, where the average ejection angle reached a maxim\am 
for 1.0 keV bombardment corresponding to the peak in 
oxygen yield. 

J. EJECTION TIME VERSUS ENERGY 

The plots of ejection time versus energy are contained in 
Appendix H. As can been seen, the highest energy atoms are 
ejected very early in the process. As time progresses, the 
ejection energy falls rapidly with a few low energy stragglers 
coming off as the collision cascade dies and recrystalliza- 
tion begins. The distributions are independent of either 
adatom density or location. The duration of the ejection 
period is again seen to be a function of the crystal atomic 
mass as was found earlier in comparing ejection time dis- 
tributions. Increasing ion energy increases the energy of 
those atoms which are ejected early, but has no appreciable 
effect on the bulk of the ejected atoms. 

The energy distribution of ejected oxygens mimics that 
of the substrate, with the most energetic atoms coming off 
early in the cascade. This close correlation is most likely 
due to fragmentation where oxygens are ejected with their 
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metal atom(s) . There iS/ though, a faster drop- 
off energy than in the case of the substrate. Some 
broadening is observed with more energetic atoms being 
ejected later as the adsorption sites are shifted to those 
involving a higher degree of bonding. 

K. SURFACE LAYER EJECTION FREQUENCY 

Appendix I contains diagrams depicting the relative 

% 

frequency of ejection of surface layer atoms. For simplicity, 
the plots contain information for a single set of 104 impact 
points. Also indicated are the locations of the impact 
zones that produced the ejections. 

The most obvious feature is the tendency for the damage 
area to extend vertically in the Z direction for the hep 
crystals, while the pattern for bcc is much more symmetric. 
This can be explained in terms of the locations of close- 
packed rows. The two close-packed rows in bcc are in the X 
and Z directions, however hep only has one of three in the 
X direction, with the other two primarily in the Z direction. 

The plots also graphically demonstrate the decrease in 
substrate yield as the oxygen density is increased, as was 
mentioned in Section B. 

L. COMPARISON WITH EXPERIMENTAL DATA 

Tables 22 and 23 provide a comparison of experimental 
data on the ratios of multimer species formed by recombina- 
tion and fragmentation with that obtained by computer 
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TABLE 22 



Comparison of Experimental and Simulation Results 
on Multimer Formation Mechanisins for Titanium 



PCPCPCPC 



RATIO 


YU 


(2x2) 

ATOP 


(2x2) 

ATOP 


(2x2) 

BR 


(2x2) 

BR 


(2x2) 

TFA+B 


(2x2) 

TFA+B 


(2x2) 

ALL 


(2x2) 

ALL 




o7t/ 


17.5 


20.0 


50.6 


21.0 


51.5 


11.7 


30.3 


15.7 


38.4 




TioVTi"*" 


50.0 


1.55 


0.60 


0.70 


0.92 


0.43 


1.43 


0.70 


1.07 


R 


TiO /Ti^ 


0.01 


7.05 


29.6 


5.40 


10.7 


1.97 


2.16 


4.00 


8.63 


F 


TiO^/Ti^ 


0.67 


- 


- 


- 


0.04 


- 


0.12 


- 


0.04 


R 


TiO”Ai^ 


0.28 


- 


- 


- 


0.36 


- 


0.18 


- 


0.26 


F 


Tio7o“ 


286. 


7.80 


1.20 


3.20 


1.80 


3.78 


4.64 


4.47 


2.78 


R 


Tio7o“ 


0.07 


35.3 


58.5 


26.3 


20.8 


17.2 


6.96 


25.3 


22.5 


F 


TiO^/O” 


3.90 


- 


- 


- 


0.08 


- 


0.44 


- 


0.09 


R 


Ti0“/0“ 


1.60 


- 


- 


- 


0.82 


- 


0.66 


- 


0.71 


F 



NOTE; All quantities are in percentages. Experimental data 
taken from Reference 22. "R" and ”F" refer to 

whether the mechanism involved was recombination or 
fragmentation. The TFA+B entries were obtained by 
averaging the data from the TFA and TFB positions. 

The ALL entires represent the average from all three 
adatom sites. 



81 



TABLE 23 



Comparison of Experimental and Simulation Results on 
Multimer Formation Mechanisms for Vanadium and Niobium 



PC PC 



RATIO 


YU 


(2x2) 

ATOP 


(2x2) 

ATOP 


RATIO 


YU 


(2x2) 

ATDP 


(2x2) 

ATOP 




0~/V^ 


10.0 


10.8 


29.1 


O’/Ito'’’ 


UNK 


13.6 


30.8 




vo'^/V’ 


27.7 


2.44 


5.19 


NbO'VNb''' 


UNK 


4.55 


5.01 


R 




0.02 


9.76 


22.1 


NbO'/Nb'*’ 


UNK 


6.04 


14.3 


F 




0.13 


- 


- 


NbO^/Nb'*’ 


UNK 


- 


0.42 


R 




0.78 


- 


- 


NbO^/Nb'*’ 


UNK 


- 


0.28 


F 


voVo” 


278. 


22.6 


17.8 


NboVo” 


UNK 


32.7 


16.6 


R 


vo“/o“ 


0.19 


90.4 


76.0 


NbO”/0” 


0.72 


43.3 


47.4 


F 


V0+/0- 


0.14 


- 


- 


NbO^o” 


UNK 


- 


1.32 


R 


VDj/O- 


7.80 


- 


- 


Nb0”/0" 


14.0 


- 


0.88 


F 



NOTE: All quantities are in percentages. Experimental data 

taken from Reference 23. ”R" and "F” refer to whether 

the mechanism involved was recombination or 
fragmentation . 
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si.mula'tion . The experimental data was obtained by Yu 
[Refs. 22/23] using 500 ev neon and argon to bombard the 
three metals. His experiments used cesiated—oxygenated 
polycrystals and through the change in the work function 
induced by the cesium, he was able to obtain information 
on the charge state dependence of the two multimer ejection 
mechanisms. The ion ratios presented in the tables were 
calculated by extrapolating the relative yield of each 
species back to a zero change in the work function which 
should represent the cesium-free case (Appendix J) . 

No exact comparison can be made due to the uncertainty 
with respect to relative ionization probability, oxygen 
placement and density, as well as differences in ion mass 
and energy. However, if the ratios of ions formed by the 
two mechanisms were similar, the computer simulations would 
support Yu's proposed model. This was not found to be the 
case. In each of the oxygen positions and densities the 
ratio of fragmentation to recombination was found to be 
opposite that reported by Yu. 

This variance between the simulation and experimental 
results is probably not due to differences in ion energy 
since Table 12 indicates that the total multimer yield 
decreases as the ion energy increases . The percentage of 
multimers formed by fragmentation also decreases, but this 
is believed due to the longer-ranged transport of energy 
which allows more ejection mechanisms to operate. The 
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activation of more mechanisms could then decrease the rela- 
tive fraction of multimers formed by fragmentation. The 
simulation yield obtained at 500 eV was, in any case, too 
small for meaningful comparison. 

Failure of the simulations to support Yu's theory may 
be caused by an inadequate treatment of the ionization proba- 
bility contribution to the determination of multimer sta- 
bility. Since the program does not determine charge state, 
some of the multimers that were classified as stable may, 
in reality, be unstable and thus reduce the percentage 
found to be formed by fragmentation. 
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V. CONCLUSIONS 



Ths program has the inherent flexiblity to permit de- 
tailed analysis of virtually any aspect of sputtering. The 
ease with which parameters may be isolated and varied 
greatly facilitates the determination of the significant 
factors involved, as well as the interdependence of these 
factors . 

The results of these simulations support earlier find- 
ings concerning the effects of channeling and that the majority 
of sputtered atoms originate in the top layers of the crystal. 
The angular distribution of ejected atoms agrees with ex- 
perimental results in the case of vanadiiam and niobium, 
however the discrepancies noted in the ejection patterns for 
titanium require further study. 

The decrease in oxygen yield from crystals with adatoms 
in the three-fold positions indicates that the oxygens are 
most stable in these sites. Such positions would also 
facilitate their incorporation into the metal lattice as the 
oxygen exposure is increased, resulting in the formation of 
a stable oxide in the top layers. It is also possible that 
the oxygens assume an A-top position at low oxygen exposures 
and as the coverage is increased, oxygen-oxygen repulsion 
causes a migration to the three-fold sites. 

The disagreement between the simulation results and 
Yu's theory and findings on multimer formation need not be 



85 



due to a failure of his basic model, but rather to an incom- 
plete understanding and treatment of the chemistry involved 
after the atoms are ejected from the surface by the simulation 
model. On the other hand it is somewhat difficult to under- 
stand how the yield of TiO^, presumably formed by recombination 
as suggested by Yu, can be so much larger than Tio”, from 
fragmentation, considering the probabilities involved in two 
or more atoms finding themselves in the same location with 
relative velocities favorable enough to permit bonding. It 
may also be that a significant percentage of the negative 
ions stabilize themselves by electron ejection prior to 
their collection in the mass spectrometer. Yu has also 
pointed out that several discrepancies in his findings may 
be due to the fact that as the oxygen density is increased 
Tio"*" and TiO^ may be formed by both mechanisms . 

An item that bears further study is the increase in the 
most probable energy of the ejected substrate atoms when 
oxygen is introduced in the A-top position. At present 
there is no apparent reason for such an increase. If this 
were a model artifact one would expect the effect to carry 
over, to some degree, to crystals with adatoms in the bridge 
positions. Instead, the average and most probcible energies 
for the bridge and three-fold sites are essentially the 
same . 
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APPENDIX A 



POTENTIAL FUNCTION GRAPHS AND PARAMETERS 
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Fig. 16. V-V solid phase potential function and associated parameters. 
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Fig. 17. Nb-Nb solid phase potential function and associated parameters 
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Fig. 18. Ti-0 solid phase potential functionand associated parameters. 
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Fig. 21. 0-0 gas phase potential function and associated parameters. 
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Fig. 25. 



97 



ATOMS PER SINGLE ION 




I g 

1 M 

ii 




ro 80 80 ro ro ro to oo 



NOUDVUd 





Fig. 26. 
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Fig. 28. 
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Fig. 30. 
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Fig. 31 . 
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Fig. 34. 
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Fig. 43. 
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Fig. 45. 
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Fig.. 47. 
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Fig. 48. 
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Fig. 50. 
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Fig. 52. 
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Fig. 53. 
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Fig. 63. 
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Fig, 64. 
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Fig.. 66. 
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Fig. 67. 
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Fig. 68. 
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Fig. 69. 
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Fig. 71. 
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Fig. 73. 



145 



200 200 
EJECTION TIME (FSEC) 





I 



f 









INITIAL LOCATIONS OF EJECTED ATOMS 




146 




147 






Fig. 76. Origin of sputtered atoms, Ti C(2x2) A-top. 
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Fig. 80 (cont.). 
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Fig. 81 (cont.). 







4. 




5. 6. 

Fig. 82. Origin of sputtered atoms, Ti C(2x2)TFB, 1.0 keV. 
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Fig. 84 (cont.). 
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Fig. 85. Origin of sputtered atoms, V and Nb clean surface. 
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Fig. 85 (cont.). Dashed lines indicate second layer atoms. 
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Fig. 86 (cont.). 
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Fig. 86 (cont.). 
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Fig. 87. Origin of sputtered atoms, V and Nb C(2x2) A- top. 
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Fig. 87 (cont.). 
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Fig. 87 (cont.). 
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Fig. 88. Ejection time difference for dimers formed from clean metal surfaces. 
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Fig. 89. Ejection time difference for Ti P(2x2) A-top. 
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Fig. 90. Ejection time difference for Ti C(2x2) A- top. 
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Fig. 92. Ejection time difference for V C(2x2) A- top. 
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Fig. 93. Ejection time difference for Ti P(2x?) bridge. 
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Fig. 94. Ejection time difference for Ti C(2x2) bridge. 
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Fig. 95. Ejection time difference for P(2x2)TFA. 
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Fig. 96. Ejection time difference for Ti C(2x2)TFA. 



1S3 





Fig. 97. Ejection time difference for Ti P(2x2)TFB. 
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Fig. 98. Ejection time difference for Ti C(2x2)TFB. 
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Fig. 99. Ejection time difference for Nb P(2x2) A-top. 
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Fig. ion. Ejection time difference for Nb C(2x2) A-top. 
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Fig, 102. 



189 






10 KEV NB(110)/AR CLEAN SURFACE 




Fig. 103. 



190 





1.0 KEV TKOOD/AR + OX P2X2TFB (TO 
SPOT PATTERN 




Fig. 104. 



1^1 



0.6 KEV *m001)/AR + OX C2X2TFB (Tl) 



SPOT PATTERN 




Fig. 105. 



192 



to KEV TKOOD/AR + OX C2X2TFB (Tl) 



SPOT PATTERN 




Fig. 106. 



19 3 



2.0 KEV TUOOD/AR + OX C2X2TFB (H) 



SPOT PATTERN 




Fig. 107. 



194 



1.0 KEV TKOOD/AR + OX P2X2 TF (A+B) 



SPOT PATTERN 




Fig. 108. 



19 5 




1.0 KEV 7U00D/AR + OX C2X2 TF (A+B) 



SPOT PATTERN 
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Fig. 112. 
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Fig. 113. 
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Fig. 114. 
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Fig. 115. 
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Fig. 116. 
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Fig^ 117. 
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Fig. 118. 
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Fig. 119. 



206 



200 too 

EJECTION TIME (FSEC) 



APPENDIX I 



SURFACE LAYER EJECTION FREQUENCY PLOTS 
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Fig. 120. Relative surface layer ejection frequency. 
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Fig, 121. Relative surface layer ejection frequency. 
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Fig. 122. Relative surface layer ejection frequency. 
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Fig. 123. Relative surface layer ejection frequency. 
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Fig. 124. Relative surface layer ejection frequency. 
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Fig. 125. Relative surface layer ejection frequency. 
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Fig. 126. Relative surface layer ejection frequency. 
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Fig. 127. (From ref. [22] ). 



214 




Work function dependence of the V*^, VO***, 
and V02'*^ yields at 5-L Oj exposure. The pnmary ion 
was 500-eV Ne"*. 




Work function dep>endence of the 0~, VO~, 
and VOi* yields. The O; exposure was 4 L, 
The pnmary ion was 5CO -«V Nc^. 




Worfe function dependence of the O", 
NbO““. Nb02*“, and NbOj* yields. The Oi exposure 
was 4 L. The primary ion was 500 -eV Ar*^. 
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